A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
01(03.13.4,12.07.1,08.16.2,08.02.3,12.04.1) 



ASTRONOMY 

AND 
ASTROPHYSICS 

February 1,2008 



A robust and efficient method for calculating the magnification of 
extended sources caused by gravitational lenses 



oo 

ON 
ON 



M. Dominik 12 * 

1 Institut fur Physik, Universitat Dortmund, D-44221 Dortmund, Germany 

2 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

Received ; accepted 



5-i , 
Oh- 

<: 

: 



> 

On 
in 
o 

o 

oo 
On 



Of 
i 1 

O 



Abstract. To determine the magnification of an extended 
source caused by gravitational lensing one has to perform a 
two-dimensional integral over point-source magnifications in 
general. Since the point-source magnification jumps to an infi- 
nite value on caustics, special care is required. For a uniformly 
bright source, it has been shown earlier that the calculation sim- 
plifies if one determines the magnification from the area of the 
images of the extended source by applying Green's theorem so 
that one ends up with a one-dimensional integration over the 
image boundaries. This approach is discussed here in detail, 
and it is shown that it can be used to yield a robust and efficient 
method also for limb-darkened sources. It is also shown that 
the centroid shift can be calculated in a similar way. 
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• i—i . 1. Introduction 

X 

5^ For fitting light curves for the ongoing microlensing events, 
?3 there is a need for robust and efficient methods for calculat- 
ing the magnification of extended sources, which are not lim- 
ited to point-lenses. Among the observed events, the presence 
of binary lenses is a reality (Dominik & Hirshfeld 1994[ 1996; 



Udalski et al. 1994; Alard et al. 1995; Bennett et al. 1996) 



and planetary events involve a special case of a binary lens. In 
addition, for some configurations, the light curve for a limb- 
darkened source will differ significantly from that of a uni- 
formly bright source. The limb-darkening effect has recently 
been observed in the galactic microlensing event MACHO 97- 
BLG-28, which involves both an extended source and a binary 
lens, by the PLANET collaboration (Albrow et al. 1998a,b); the 
fitting has been done by myself using the algorithm described 
in this letter. 

If one wants to integrate the point-source magnification in 
two dimensions one has to take special care of the position 
of the caustics, where the point-source magnification becomes 
infinite. While this integration can be performed easily for a 



1994; Sahu 1994; Dominik 1996), this would be a dif 



1992, p. 313; Witt & 



ficult task for a general lens (e.g. a binary lens), especially at 
a cusp singularity. In contrast, the area of the images of the 
extended source and therefore its magnification remains con- 
tinuous when the source hits a caustic. The determination of 
the extended source magnification from the boundaries of the 
image areas has been used to analyze the image s of ba ckground 
galaxies behind a cluster of galaxies (Dominik 1993 ). The im- 
age boundaries can be obtained with a contour plot of an im- 
plicit function describ ing th e source boundary in the lens plane 
(Schramm & Kayser 1987 ). This method has been expanded 
with routines for correcting, testing and finally analyzing the 
contour line in order to produce an efficient and safe algorithm 
(Dominik 1995). In that paper, it is noted that it is easy to an- 
alyze the images from the contour line data, and an example 
is given, where quantities such as the area, width, length and 
curvature of the image have been determined. 

Concerning microlensing light curves, it has been noted by 
Bennett & Rhie ( 1 996 ) that it is advantageous to integrate in the 
lens plane rather than in the source plane to determine the mag- 
nification of an extended source. For uniformly bright sources, 
Gould & Gaucherel (1997) proposed applying Green's theorem 
so that only one integration along the image boundary must be 
performed rather than two over the image area. This approach 



is identical to that used earlier (Dominik 1993 1995). The con 
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tour plot method is the most convenient way to obtain data 
points on the image boundary from which the area can be cal- 
culated. In Sect. 2, this general approach is described, Sect. 3 
gives details for a uniformly bright source, and Sect. 4 shows 
how this approach can also be used for limb-darkened sources, 
in which case an easy-to-perform two-dimensional integration 
remains. In Sect. 5, the calculation of the centroid shift is dis- 
cussed. 



2. Magnification and Green's theorem 

The magnification is given by the ratio of the area of the images 
to the area of the source in the absence of the lens. Let (yi, j/2) 
denote arbitrary cartesian source coordinates and (0:1,352) de- 
note corresponding image coordinates, where in the absence 
of the lens y = x. Let I s (y) denote the surface brightness of 
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the source and h(x) denote the surface brightness of the im- 
ages. The conservation of photon number then requires (see 
e.g. Schneider et al |1992| , p. 33; Schramm & Kayser |1987 ; 
Kayser & Schramm |1988[) that 



h{x) = I B (y(x)), 



(1) 



where y(x) gives the source position y related to the image 
position x. 

With A s being the region in the (yi, j/2)-plane subtended 
by the source and A\ being the region in the (x\, X2)-plane 
subtended by the images, the magnification is given by 



j Ai h{x)dxi dx 2 



(2) 



In the limit of a point source, one obtains \i = J] P>(xj), 
where the sum runs over the m images Xj of the source, and fi 
is given by 



Jx{Xj) = 



1 



det 



(£)(« 



(3) 



The normalization of the brightness profile function I s can 
be chosen so that the integral in the denominator (Eq. (Q)) be- 
comes the area of the source, and for a circular source of radius 



4(y) dyi dy 2 = 7ri? s 2 



For a uniformly bright source, one has I s (y) = 1 for posi- 
tions within the source, so that h(x) — 1 for positions within 
the images, and the magnification becomes 



ttR 2 ,., 



dxi dx2 



(5) 



where the remaining integral is just the area of the images. 

Green's theorem now states that for two functions 
P{x-y,X2) and Q{x\,X2) which are continous in a region A 
and whose partial derivatives dP/dx2 and dQ/dxi are also 
continuous in A, 



dQ dP\ 

o a — dxi dx 2 

ax\ 0x2 J 



P dx\ + Q dx2 



(6) 



OA 



where dA denotes the boundary of the region A, which has to 
be piecewise-smooth. Note that the values of P and Q outside 
the region A do not play a role and can be chosen arbitrarily as 
long as P and Q satisfy the conditions above. 

3. Uniformly bright sources 



For a uniformly bright source, one can now choose P = 
and Q = ^x\ to obtain 



1 



2-KRt. 



X\ dx2 — X2 dx\ 



\X2 



(7) 



Fig. 1. Two successive points a;W and x^ l+1 ^ on the image 
boundary and the corresponding area segment 

For n discrete positions ajW on the image boundary, where 
x (n+i) _ x (i) t w hich are obtained e.g. by a contour plot, the 
magnification can be approximated by 



1 



'■ IC . 1 



4 l) ) 



(8) 



This means that one adds up rectangular segments with the 



length x± and the width Ax^ 



„(*+!) 



(and corre- 

(i+l) 



sponding segments by interchanging the axes). However x\ 
has an equal "right" to be convolved with Ax2(i), yielding the 
symmetric version 



(4) n = 



1 



i=l 



(4 ] 



4 i+1) )(4 i+1) 



(9) 



which sums over trapezoidal segments having an area which is 
the mean of the rectangular areas formed from segments with 
length Xi and Xj* + , and corresponds to a replacement of the 
true boundary by a polygon with n corners at x^ (see also 
Figure [j]). By increasing the number of points used, the magni- 
fication can be determined to the desired precision. 

4. Limb-darkened sources 

Let us now consider a limb-darkened source at y^ with radius 
i? src , with the brightness profile 



Uy) = h 

where 

fi(r;u) 



(y-y (0) ) 

R 2 



1 



1 



l-u/3 

for < r < 1, and fi(r,u) = for r > 1. This profile 
is normalized so that Eq. (Q) is satisfied. The parameter u is 
chosen from the interval [0, 1], where u = corresponds to a 
uniformly bright source. This brightness profile is the simple 
'linear' modelQ widely used, though for some type of stars, one 

1 The profile involves a linear term of cos 6, where 9 is the angle 
between the normal to the surface and the direction to the observer. 



(10) 



(11) 
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should add other terms (e.g. Claret et al. 1995). The general 
ideas of the approach discussed here do not depend on this spe- 
cial choice, especially terms like (1 — r 2 )™ (i.e. like cos 2 ™ 8) 
can be treated in complete analogy to the \fl — r 2 -term.^| 
From Eqs. O) and (0), one obtains the magnification as 



M = n 9 / h{x)dxi dx 2 ■ 

*Rirc J A, 



(12) 



To apply Green's theorem one has to find functions P and 
Q which satisfy 



dQ OP 

a a — = Il{Xi,X 2 ) = I s {y(Xi,X 2 )) . 

ax% 0x2 

Such functions are given by 

I r x 2 



P(xi,x 2 ) 
Q{xx,x 2 ) 



2 ./.. 

■Xl 



Is {y{xi,x' 2 )) dx' 2 , 
I s (y{x[,x 2 )) dx[ . 



(13) 

(14) 
(15) 



Since the brightness profile I s and the lens equation y(x) are 
given in analytical form, the integrals can be evaluated numer- 
ically in general. 

Since the brightness profile has an infinite slope at the limb 
of the source, it is advantageous to use another continuation for 
r > 1.0 Let us write the brightness profile in the form 



h{r; u) = C{u) (uB(r) + l-2u) , 
where 

C(u) ' 
and 



l-u/3 



B(r) 



for < r < 1 
for r = 1 



(16) 
(17) 

(18) 



1 - y/l - 1/r 2 for r > 1 



This definition is identical to the previous one except for r > 
The functions P and Q read 



P(xi,x 2 ) 



-C/2 



{l-2u)(x 2 -x 2 0) ) + 



+ u B(r(xi, x' 2 )) dx' 2 



(19) 



Q(x 1 ,x 2 ) = C/2 



(l-2u) (.ti-x^ 0) ) + 



+ u 



(0) 



B(r(x' 1 ,x 2 ))dx' 1 



(20) 



2 In fact, for the discus sion of the MACHO 97-BLG-28 event (Al- 
brow et al. 1998a,b), a v^l — r 2 -term has also been included in the 
brightness profile. 

3 As mentioned earlier, the value for r > 1 does not play a role. 

4 One may also use other continuations for r > 1 which decrease 
differently for r — > 00. 



where 



r(xx,x 2 ) 



(y{xi,x 2 ) - yW)' 
R 2 



(21) 



The function B(r) has been chosen, so that for all r, B{r) > 
in order to avoid contributions of different sign in the numerical 
integration process. The tail of the function is limited to 



/>oo 

/ B(r) dr = ^ - 1 w 0.571 , 
Ji 2 



(22) 



so that the integral is not dominated by the tail contribution. 

One remaining point of interest is the choice of the lower 
integration bound a;^ ' . To avoid unnecessary integration in the 
tail region of B(r) that would result for a fixed choice for 
such as (0, 0), the lower integration bound can be chosen as the 
center of each image which can be determined as shown in the 
next section. 

With these funtions P and Q one can approximate the mag- 
nification by the expression 



1 - 
OttR 2 



Q ( x l ' 2 ^2 ' 



) + 



„(') , Ji+V- 



(x. 



p[^+x^),4 +1) 



(23) 



which uses the symmetries and reduces to Eq. (^) with = 
(0, 0) for the case u — (uniformly bright source). An example 
for a light curve for a limb-darkened source and one with a 
uniformly bright source behind a binary lens is shown in Fig. 2. 

5. The centroid shift 

The position of the centroid of light can be determined in a 
similar way. For m images of the source, the position of the 
centroid of light x c is given by 



J A . xl\(x) dxi dx 2 
3=1 J 

m 

J2 J A . h(x) dxi dx 2 



(24) 



where an integration has to be performed over the regions sub- 
tended by the images Aj . The integral in the denominator is just 
that discussed in the previous section and the integral in the nu- 
merator differs only by an additional factor of x. The integral 
can be solved in the way described in the last section by only 
replacing I\(x) by Xi I\(x) (i — 1, 2). For a uniformly bright 
source, the corresponding functions P and Q are 



P(xi,x 2 ) 



1 

--XIX2 



Q(xi,x 2 ) = \x\ 



(25) 
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Fig. 2. The light magnification for a 
limb-darkened source with u = 1 
(curve with larger peak amplification) 
and a uniformly bright source behind a 
binary lens. The angular source radius 
corresponds to 0.03 #e, where 6*e de- 
notes the angular Einstein radius 5 . The 
mass ratio between the two lens com- 
ponents is 4, their angular separation 
is 0.68 #e- The source center passes 
perpendicularly to the line connecting 
the lens objects at a minimal impact 
of 0.004 6*e from the midpoint towards 
the lighter component. £e is the time in 
which the angular source position rela- 
tive to the lens changes by 9e, t = 
corresponds to the point of time when 
the center of the source crossing the 
line connecting the lens components. 
The inset shows the caustics, the posi- 
tion of the lens objects (crosses), and 
the size and trajectory of the source 
which sweeps over a cusp caustic; the 
coordinates are multiples of 9e- 



for the a; i -component and 
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P(x 1 ,x 2 ) = , Q(xi,x 2 ) = -Xix 2 

for the X2 -component, so that one obtains with Eq. 
expressions 



(26) 
the 



x c ,i 



(X 



(i+1) 



+ ( x [ i) -x^ 1 ' ){x 



„(*+!) 



>}■ 



1 

^ 2 = -g I E{( 

2 — 1 

+ (4 l)a +4 w,a )(*r , -*? ) )}, 



4 i)2 -4 i+1)2 )(4 i+1) +4 i) ) + 

i2. 



(27) 
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where A denotes the area of the image, which can be deter- 
mined as shown in Sect. [|. 

For a point source, the centroid's position can be written as 



(29) 



If one compares this expression with that given by Eq. ( 24 ), one 
sees that the (finite) source size has cancelled out in Eq. (24j). 
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For a lens at a distance and a source at a distance D s from 
the observer, the angular Einstein radius for a lens of total mass M is 

given by = (4GMc~ 2 D" 1 D' 1 (D s - 7J d )) 1/2 . 
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Figure 1 
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Figure 2 
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